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Abstract 

We review our analysis of the strong coupling limit of compact QED on a 
lattice with staggered Fermions. We show that, for infinite coupling, compact 
QED is exactly mapped in a quantum antiferromagnet. We discuss some 
aspects of this correspondence relevant for effective field theories of Josephson 
junctions arrays. 
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1 Introduction 



The strong coupling limit of gauge theories such as QCD exhibits many of the 
qualitative features of their spectrum. For example, confinement is manifest 
and also, in some cases, it is straightforward to show that chiral symmetry 
is broken. Recently this issue has been re-examined in detail [|I], and the 
possibility for chiral symmetry breaking for a wide variety of gauge theories 
has been explored. 

It has been recognized for some time that the strong coupling limit of 
lattice gauge theory with dynamical fermions is related to certain quantum 
spin systems. This is particularly true in the Hamiltonian picture and ap- 
peared in some of the earliest analyses of chiral symmetry breaking in the 
strong coupling limit [^. In that work, the chiral symmetry breaking phase 
of strong coupling gauge theory was analyzed by drawing an analogy with 
quantum antiferromagnets and then using mean field theory and spin-wave 
analysis to compute the properties of the pion spectrum. 

In more recent literature it has been noted that there are several for- 
mal similarities between some condensed matter systems with lattice fermions 
(particularly certain antiferromagnetic spin systems) and lattice gauge the- 
ory systems, usually in their strong coupling limit. For example, it is well 
known that the quantum spin 1/2 Heisenberg antiferromagnet is equivalent 
to the strong coupling limit of a particular U(l) lattice gauge theory [Q. It 
can also be written as a kind of SU(2) lattice gauge theory P]. 

An important feature of the relationship of lattice gauge theory with spin 
systems is the fact that the staggered Fermions which are used to put the 
Dirac operator on a lattice resemble ordinary lattice fermions used in tight 
binding models in condensed matter physics when the latter have a half-filled 
band and are placed in a background U(l) magnetic field vr (mod 27r) (1/2 
of a flux quantum) through every plaquette of the lattice. This is an old 
result for d = 2 + 1 where it was already recognized in the first work on 
the Azbel-Wannier-Hofstaeder problem and has since been discussed in the 
context of the so-called flux phases of the Hubbard and quantum Heisenberg 
models. It is actually true for all d > 2 + 1 and this equivalence was used in 
[0, 1^ to find mappings between strong coupling lattice gauge theories with a 
wide variety of fermion contents and various quantum spin systems. 

A place in condensed matter where the magnetic flux which is necessary to 
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produce a relativistic spectrum for the fermions can be achieved is in the flux 
phases 0. In this case the magnetic field comes from a certain condensate. 

Alternatively, as an external magnetic field, for ordinary lattice spacings, 
it is as yet an experimentally inaccessible flux density. Flux densities of this 
order of magnitude can presently only be achieved in analog experiments 
where macroscopic arrays of Josephson junctions take the place of atoms 
at lattice sites. These systems are described in mean field theory by tight 
binding electrons in a background U(l) magnetic field P, |^ with accuracy 



increasing with d For d = 2 their ground state is expected to be a flux 
phase [|1^, |I2[ . The low energy limit in this phase resembles strongly coupled 
QED. 

In the following we shall review some of the features of strongly coupled 
QED which are useful in mapping the strong coupling limit onto spin systems. 
We shall also emphasize some of those features which we feel are valuable 
in formulating effective field theories for Josephson junction arrays. Let us 
begin with a brief summary of the physics of the latter arrays. In the classical 
limit they are described by the effective Hamiltonian 

H = -J cos{e, - e, + ai,) (1) 

where 6i are phases of the superconducting order parameter at the i'th lat- 
tice site (which is an island of superconducting material joined by a weak 
connection to other similar islands), atj is proportional to the line integral 
of the electromagnetic vector potential along the connection between sites i 
and j, 

dij = — A - al 
he Ji 

Also, J is a coupling constant which is generally a function of external mag- 
netic field and temperature and the summation is over nearest neighbors. 

The uniformly frustrated system occurs when the external magnetic field 
is a constant in each plaquette (where we assume that the islands lie on a 
two dimensional lattice), 

E«^. = 2vr/ (2) 

□ 

and / is the ratio between the magnetic flux threading the plaquette and the 
unit flux. The latter system can be mapped onto a classical Coulomb gas 
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problem with Hamiltonian |n 

H = -7rJ^[m(r) + f]G{r - r')[m{r') + f] (3) 

r,r' 

where the charges m(r) reside on the dual lattice, are restricted to integer 
values (m G 2) and the total charge is constrained as 

- — - m(r) = f 

{\V\ is the total number of lattice sites). In D = 2 the charges represent 
vortices in the phase field of the superconducting order parameter. Even in 
that case, the ground state of this system is known only for some particular 
filling fractions / which depend on the lattice. For example, on a square 
lattice the ground states are known for some rational filling fraction / = |, 



p, g G 2 |]T0|, Q . Of particular interest to us is the fully frustrated case, f = \ 



on a square lattice. There, the ground state is a staggered configuration with 
alternating m + f = \ and m + / = — ^ on neighboring sites. 

Also, for this charged plasma in the continuum, it is known that the 
ground state of a classical Coulomb gas of this kind is a triangular Wigner 
lattice [l^. Furthermore, when f = \ there is a U(l) symmetry and the 



ground state exhibits a discrete symmetry 7L2. At zero temperature the 
system chooses one of the two phases. 

It has also been shown that the equation which determines the configu- 
ration of the order parameter is equivalent to the equation which determines 
the spectrum of electrons in a lattice in a background magnetic field, the 
so-called Azbel-Wannier-Hofstaeder problem p, P). 

These classical systems are the limit of a quantum system where the 
variables conjugate to the phases of the order parameters are the electric 
charges residing on the islands 

[9,, Qj] = 2ieh5ij (4) 

The effective Hamiltonian is 

H =]: Yl{Qi - Qxi)C~j^{Qj - Q^j) - J cos(0i - 4>j + ttij) (5) 

ij <i,i> 
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where Qi is the electric charge and Q^i is an offset, or baclcground charge 
residing on island i. The Coulomb interaction of the charges is described 
by the capacitance matrix Cij. The classical Hamiltonian is regained in the 
limit Cij — s> oo or J — s> oo. 

Another model which is often compared with @, and in fact which has 
been argued to describe the critical behavior of non-classical junction arrays 
[ |T5[] , is the Bose-Hubbard model with Hamiltonian 

H = -'^ E (<f.V'^-$,+h.c.) + i5:i-,,n,(n,-5,,)-/z^n, (6) 

<i,j> ij i 

where $j and $| are the annihilation and creation operator for a boson at 
site i, respectively and rii = is the boson number at site i. The total 
number of bosons is controlled by the chemical potential /i < 0. The on-site 
part of the self-interaction vanishes if the site is occupied by a single boson. 
The energy for more than one boson on a site should be large, to give the 
bosons a hard core. With hard core bosons, this model is equivalent to the 



frustrated quantum XY model [12 



There is also a mapping of the frustrated quantum XY model onto a 
lattice Chern-Simons theory coupled to Fermions [|T^]. For the uniformly 
frustrated case, the mean field theory of the latter system is equivalent to 
lattice fermions in an external field, related to the frustration parameter. 
The fermion spectrum is that of the Hofstaeder problem and is similar to the 
fermion system encountered in the flux phase of the Hubbard model. 

It is an interesting question to what extent these models share other fea- 
tures of a lattice field theory such as quantum electrodynamics. If the fermion 
spectrum is relativistic and, given that the models have a U(l) gauge invari- 
ance, one might expect dynamical effects to generate the relevant operators 
of electrodynamics. Furthermore, the infinite coupling limit of lattice QED 
reduces to the problem of solving for the ground state of a classical coulomb 
gas with a logarithmic potential, identical to the classical limit of the junc- 
tion array problem. The ground state of this coulomb gas system has been 
examined in both contexts in previous literature. 

It is interesting to speculate that, when we move away from the classical 
limit in the case of junction arrays, and when we move away from the infi- 
nite coupling limit in the case of compact QED, we still have systems with 
many common features. For example, both systems have a first order phase 
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transition. In the case of the juction array it is the transition which destroys 
superconductivity. In the case of compact electrodynamics it is where chiral 
symmetry^, which is broken in the strong coupling limit, is restored. It is 
accompanied by a confinement-deconfinement transition for electric charges. 
In the coulomb gas representation of the classical junction array the phase 
transition corresponds to melting of the Wigner lattice to form a two di- 
mensional plasma. In compact electrodynamics the phase transition restores 
full translation invariance of the lattice theory (corresponding to chiral in- 
variance in the continuum theory) and the electrons are deconfined. In the 
latter case the strong couphng, confining ground state has the structure of a 
Wigner lattice of electrons. 

In the following we shall review some aspects of strongly coupled lat- 
tice QED relevant to an effective field theory of Josephson junction arrays. 
We shall show how the analogy between Dirac fermions and tight-binding 
condensed matter fermions is used to find an exact mapping of the strong 
coupling limit of compact lattice QED with dynamical staggered fermions 
onto certain quantum SU(A^) antiferromagnets. The precise structure of the 
resulting antiferromagnet depends on the number of fermion flavors. We 
shall also review the argument by which QED with staggered Fermions 
in any number of spatial dimensions, D > 1, resembles a condensed matter 
system of lattice electrons in an external magnetic field and with a half-filled 
band . 

Of course, the understanding of the phase structure, particularly the con- 
fining and chiral symmetry breaking phases, and nature of phase transitions 
of QED is of great fundamental interest to particle physicists. These issues 
have been studied in the continuum quantum electrodynamics where the 
possibility of a phase transition in the strong coupling regime is suggested 
by approximate analytical results using the quenched ladder approximation. 
The idea for the mechanism, which has long been advocated by Miranski and 
collaborators [ITB, [T^, is that when the electric charge reaches a certain criti- 
cal value the vacuum of QED becomes unstable, leading to a sort of collapse 
phenomenon, analogous to the sparking of the vacuum in the presence of a 
supercritical nucleus. In the picture which we propose, the analog of the col- 
lapse phenomenon is the formation of either charge or isospin density waves 



^We shall see in following sections that, for staggered fermions, chiral symmetry corre- 
sponds to invariance under translations by one lattice site 
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and the resulting reduction of the lattice translation symmetry from trans- 
lations by one site to translations by two sites. It occurs when the exchange 
interaction of the electrons, which is attractive, dominates the tendency of 
the direct Coulomb interaction and the kinetic energy to delocalize charge, 
giving an instability to the formation of commensurate charge density waves. 
This forms a gap in the Fermion spectrum and a particular mass operator 
obtains a vacuum expectation value. This gives an intuitive picture of how 
strong attractive interactions in a field theory can form a coherent struc- 
ture. Here, the commensurate density waves in the condensed matter system 
correspond to a modulation of the vacuum charge or isospin density at the 
ultraviolet cutoff wavelength in the field theory. 

In the continuum, we are generally interested in massless quantum elec- 
trodynamics with the action 

S = J d^+^x i-^F,,F'^^ + E + A,)r\ (7) 

where A is the ultraviolet cutoff, e is the dimensionless electric charge and 
there are Np flavors of 2'*^ ''^'-dimensional Dirac spinors. (Here, [(-D+l)/2] 
is the largest integer less than or equal to {D + l)/2.) In even dimensions 
(when D+1 is even), the flavor symmetry is S\Jfi{Np)xS\JL{Np). In odd 
dimensions there is no chirality and the flavor symmetry is SI]{Nf). What is 
usually referred to as chiral symmetry there is a combination of parity and a 
Z2 subgroup of the flavor group. We shall use a lattice regularization of (|^) 
and study the limit as 00. We shall use Np flavors of staggered fermions. 
In the naive continuum limit, this gives Np = Nl 2-component fermions 
in D = 1, Np = 2Nl 2-component fermions in D = 2 and Np = 2Nl 4- 
component fermions in D = 3. Though the lattice theory reduces to (|^) in the 
naive continuum limit, the lattice regularization breaks part of the flavor sym- 
metry. All of the operators which are not symmetric are irrelevant and vanish 
as the lattice spacing is taken to zero. In D = 3, the S\J ji{Np)xS\J i{Np) 
symmetry is reduced to S'[J{Np/2) (in D = 3, Np = Np/2) and translation 
by one site in the 3 lattice directions. These discrete transformations corre- 
spond to discrete chiral transformations in the continuum theory. \n D = 2, 
the lattice has SU(A^i?/2) (in D = 2, Np = Np/2) symmetry and two discrete 
("chiral") symmetries. In D = 1, the Sl]L{Np)xSl]p{Np) symmetry is re- 
duced to S\J{Np) (in D = 1, Np = Np) and a discrete chiral transformation. 
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In each case, the discrete chiral symmetry is enough to forbid Fermion mass 
and it is the spontaneous breaking of this symmetry which we shall examine 
and which we shall call "chiral symmetry breaking" . We shall also consider 
the possibility of spontaneous breaking of the SU(A'^p) flavor symmetry. In 
the continuum the SU(A^i;') corresponds to a vector-like symmetry subgroup 
of the full flavor group. Since the lattice and continuum theories have differ- 
ent symmetries, the spectrum of Goldstone bosons, etc. would be different 
in the two cases. We shall not address the source of these differences but 
will define QED by its lattice regularization and discuss the realization of 
the symmetries of that theory only. 

In the strong coupling limit, lattice QED with A^l lattice flavors of stag- 
gered Fermions is exactly equivalent to an SU(A^l) quantum antiferro- 
magnet where the spins are in a particular fundemental representation of 
the SU(A^l) Lie algebra. Furthermore, mass operators of QED correspond to 
staggered charge and isospin density operators in the antiferromagnet. Thus, 
the formation of charge- density waves corresponds to chiral symmetry break- 
ing and the dynamical generation of an iso-scalar Fermion mass whereas Neel 
order of the antiferromagnet corresponds to dynamical generation of Fermion 
mass with an iso-vector condensate. As a result of this correspondence we 
are able to use some of the properties of the quantum antiferromagnets to 
deduce features of strongly coupled QED. 

In the infinite coupling limit, the properties of the electronic ground state 
of compact and non-compact QED are identical. The fact that compact QED 
confines and non-compact QED (at least in high enough dimensions) does 
not confine affects only the properties of the gauge field wavefunctions and 
the elementary excitations. In the following we shall focus our attention only 
on compact QED. 

An interesting result of our approach [0] is that the case of even A^l and odd 
A^L are very different: 

When A^L is odd, the vacuum energy in the strong coupling limit is pro- 
protional to e^, the square of the electromagnetic coupling constant. We 
also find that chiral symmetry is broken in the strong coupling limit for all 
odd A^L > 1 and for all spacetime dimensions D + 1 > 2. The mass operator 
which obtains a nonzero expectation value is a Lorentz and flavor Lie algebra 
scalar. There are also mass operators which are Lorentz scalars and which 
transform nontrivially under the flavor group which could get an expectation 
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value and break the flavor symmetry spontaneously if A'l is small enough. 
When TVl increases to some critical value there is a phase transition to a 
disordered phase WE . 



In contrast, when A^l is even we find that the vacuum energy in the 
strong coupling limit is of order 1. We find that chiral symmetry is broken 
when the spacetime dimension D + 1 > 3 and when N-^ is small enough. 
The mass operator which gets an expectation value is a Lorentz scalar and 
transforms nontrivially under the flavor SU(A^l)- Thus, flavor symmetry is 
spontaneously broken. As in the case of odd A^l, there is an upper critical 



A'l where there is a transition to a disordered phase . 

For A^L odd finding the gauge invariant electronic ground state of strongly 
coupled lattice QED is equivalent to finding the exact ground state of the 
generalized classical Coulomb gas model in D-dimensions with Hamiltonian 

^coui = ^ II P{x)g{x - y)p{y) (8) 

<a:,j/> 

where the variable p{x) lives at the sites of a cubic lattice with spacing one 
(with coordinates (xi, . . . , x/j) and integers) and takes on values 

- — Y ' ■ ■ ■ ' ^ ' ^ A/l an odd mteger (9) 

and the interaction is long-ranged 

g{x — y) \x — y\'^^^ , as |x — ?/| — > oo (10) 

In any space dimension, D > 1 we prove that there are two degenerate ground 
states which have the Wigner lattice configurations 

p(x) = ±^(-i)i:"i^'' (11) 

When A^L = 1 this is a long-ranged Ising model which is known to have 
antiferromagnetic order, even in one dimension. It is also in agreement with 
a known result about the Ising model in two dimensions [0 . 

The Coulomb gas representation of the classical Josephson junction array 
is related to this system when / = 1/2, i.e. in the fully frustrated case. There, 
the charge on a site is only limited by dynamics, whereas in electrodynamics it 
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is restricted by the Pauli principle to be between —Nl/2 and Nl/2. However, 
in the case of the classical Coulomb gas, the Coulomb energy of more than 
one charge residing on a lattice site is infinite, so effectively, whenever / < 1 
the charges at lattice sites should also be either 1,0 or —1. Therefore the 
problem of determining the charge (vortex) distribution in a fully frustrated 
junction array is equivalent to finding the ground state of infinite coupling 
QED with Nl = 1. It will turn out that the solution of this problem is 
similar for all odd A'^^- 

The difference between even and odd A'l can be seen to arise from a 
certain frustration encountered whenQ, with odd A^l, one simultaneously im- 
poses the conditions on lattice QED which should lead to charge conjugation 
invariance, translation invariance and Lorentz invariance of the continuum 
limit. This frustration is absent when A^l is even. This is not an anomaly in 
the conventional sense of the axial anomaly or a discrete anomaly encountered 
in the quantization of gauge theories, as no symmetries are incompatible, but 
it is nevertheless an interesting analog of the anomaly phenomenon. 

Note that in 2+1 dimensions the difference between even and odd A^l is 
the difference between an even and odd number of 4-component Fermions. 
Our result that for an odd number of 4-component Fermions the chiral sym- 
metry is always broken for large coupling seems to contradict the continuum 
analysis in |19|, We do not fully understand the reason for this, but specu- 
late that it is related to the lattice regularization. Their model is very similar 
to ours in that they effectively work in the strong coupling limit when they 
replace the ultraviolet regularization which comes from having a Maxwell 
term in the QED action by a large momentum cutoff. In our case we have a 
lattice cut-off and the result should be very similar. Note that we agree with 
them when A^l is even and there is an even number of 4-component Fermions. 
In that case, we are also in qualitative agreement with recent numerical work 
[21] which, since it uses Euclidean staggered Fermions and there is a further 
Fermion doubling due to discretization of time, can only consider the case 
where there is an even number of 4-component Fermions. The anomalous 
behavior that we find with odd numbers of 4-component Fermions exposes a 
difficulty with treating the flavor number A^ as a continuous parameter. 

A hint as to why 2+1-dimensional Fermions should come in four-component 



^Note that this frustration is different from that caused by the external magnetic field 
in the junction array 
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units appears when we formulate compact QED in the continuum using the 
SO (3) Georgi-Glashow model with spontaneous breaking of the global sym- 
metry, S0(3)^U(1) and the limit of large Higgs mass We begin with 
the model 



The spectrum contains a massless photon which in the A — >■ oo and t;^ — >■ oo 
is the only light excitation. Since the U(l) gauge group is a subgroup of 
S0(3), it is compact. We wish to make 2+1-dimensional electrodynamics by 
coupling this model to Fermions. It is known that, if the resulting theory 
is to preserve parity and gauge invariance simultaneously, we must use an 
even number of two-component Fermions which, in the minimal case are also 



Here, the Fermions are four-component and have a parity invariant mass and 
Higgs coupling. Each SU(2) doublet contains two electrons (which can be 
defined so that both components of the doublet have the same sign of electric 
charge), thus the basic Fermion for compact QED has eight components. The 
maximal chiral symmetry is obtained in the massless case with vanishing 
Higgs coupling. It is possible, by suitable choice of mass and Higgs coupling 
i.e. m = ±g <</)>, (and also reduction of the chiral symmetry) to make four 
of the Fermions heavy, leaving four massless components. This is consistent 
with parity and gauge invariance. It is also interesting to note that staggered 
Fermions on a Euclidean lattice, where time is also discretized, produce eight 
component continuum Fermions. 

In Section 2 we review the essential features of the staggered fermion 
formalism. In Section 3 we elucidate how spin algebras may be represented 
by Fermion bilinears. In Section 4 we formulate lattice QED and analyze 
its discrete symmetries. Section 5 is devoted to the analysis of the strong 
coupling limit of compact QED. There we shall review the properties of the 
electronic ground state for both TVl even and odd. Section 6 is devoted to 
concluding remarks. 



^ = jr2 E(^;.)' + 7^ ■ wr + t((E w - ^"f (12) 

a=l a=l ^ a=l 



SU(2) doublets [|| 




(13) 
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2 Staggered Fermions 



Here we make explicit the fact that latticized (staggered) relativistic fermions 
are equivalent to the Hofstaeder problem with flux density 1/2 flux quanta 
per plaquette. 

In the following, we use units in which the lattice spacing, the speed of 
light and Planck's constant are all equal to one. Space is a D-dimensionsal 
cubic lattice with lattice points x = {xi,X2, ■ ■ ■ ,xd) {xi integers), oriented 
links [x, i] labeled by their site of origin and their direction, i = ±1, . . . , ±D, 
and oriented plaquettes [x, with corners x,x + i,x + i + j,x + j and 
orientation i x j where i is the unit vector in i-direction. 

In relativistic field theory, the purpose of the staggered Fermion method 
is to minimize Fermion doubling which always occurs when one tries to put 
Fermions on the lattice. Generally, staggered Fermions are obtained by the 
spin-diagonalization method. To implement this method, we begin with the 
naively latticized Dirac Hamiltonian for free relativistic fermions, 

Hf = ^Y1 {i^Hx)ia^Vjij{x) - {Vjij^){x)ia^ij{x)) (14) 

where the lattice derivative operator is defined as difference operator 

VJ{x) = f{x + l)-f{x) 

and are the 2[*^'^+^-'/^l-dimensional Dirac a-matrices. (Here [{D + l)/2] is 
the integer part of (-D + l)/2.) The latter are Hermitean, a^'^ = a\ and obey 
the Clifford algebra 

{a\a^} = 26'^ (15) 

They are therefore unitary matrices, a^^a* = 1. Using the properties of the 
difference operator, (|1^) can be presented in the form 

Since the Dirac matrices are unitary, the naive lattice Fermion Hamiltonian 
in (|16|) resembles a condensed matter Fermion hopping problem with a back- 
ground U(2[^^+^)/^]) gauge field given by the a-matrices. In any plaquette of 
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the lattice, [x,i,j], this background field has curvature 

a'a^a'^a^^ = -1 (17) 

The curvature resides in the U(l) subgroup of 11(2'*^ '^^^-'/^l) and has exactly 
half of a U(l) fiux quantum per plaquette. This is true in any dimensions. 
We observe that either 1/2 or zero fiux quanta are the only ones allowed 
by translation invariance and parity and time reversal symmetries of the 
Hamiltonian. 

Since the curvature of the a-matrices is U(l)-valued, we should be able 
to do a gauge transformation which presents the matrices themselves as U(l) 
valued gauge fields (i.e. diagonal). A specific example of such a gauge trans- 
formation due to Kluberg-Stern et.al. 0| is 

ij{x) {a^' {a^^ . . . (a^y^^ix) (18) 

Then 

ij\x + j)a^ilj{x) {-l)^i'=i''^-^^{x + 3)ilj{x) (19) 
The resulting Hamiltonian is 

Hp = -^E(-l)^'*^'^"' {iJ^{x+j)ij{x)-ij\x)i:{x + j)) (20) 

This describes 2[(^+i)/2] identical copies of Fermions with the same Hamil- 
tonian which must all give Fermions with the same spectrum as the original 
Hamiltonian in (|14|). Staggered Fermions are obtained by choosing one of 
these copies. This reduces the Fermion doubling by a factor of the dimension 
of the Dirac matrices, 2^^^'^^^^'^^ 

In the staggered Fermion method, we treat the components of the original 
lattice Dirac Hamiltonian as fiavors, rather than components of the relativis- 
tic spinor necessary for Lorentz invariance. The spinor components now 
reside on adjacent lattice sites. In this method, the continuous chiral sym- 
metry of the massless Hamiltonian, under the transformation ip e*^ is 
lost. There is a discrete chiral symmetry, corresponding to translations by 
one lattice site in any direction. Explicitly, 

4j{x) +j) (21) 
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is a symmetry of the Hamiltonian (|20| ) and corresponds to a discrete chiral 
transformation. 

Mass operators correspond to staggered charge densities. The operator 
S = ^(-1)^^=1 (22) 



changes sign under the chiral transformations (^T]) and corresponds to a 
certain Dirac mass. 

With staggered Fermions there is still a certain amount of Fermion dou- 
bling. The doubling can be counted by noting that the staggered Fermion 
Hamiltonian (pOf ) is invariant under translations by two lattice sites. There- 
fore, a unit cell is a unit hypercube of the lattice, containing 2^ sites and 
staggered Fermions correspond to a 2^ component spinor. The dimension 
of the Dirac matrices is d^^^'^^^^'^^. Therefore the number of Dirac spinors 
we obtain is 2^/2^^^^^^^"^^. For lower dimensions the minimum number of 
continuum flavors can be tabulated as 

d dim. of Dirac matrices No. of flavors 

1 + 1 2 1 

2 + 1 2 2 

3 + 1 4 2 

Only in 1+1- dimensions do we get a single species of Dirac Fermion. 

To see how to take the continuum limit explicitly, consider the case of 
d = 3 + 1. There, we divide the lattice into eight sublattices and label the 
spinor components as 

(23) 
(24) 
(25) 
(26) 

In terms of these spinors, the Hamiltonian (^) can be written as the matrix 
operator 

H= f d^k^\k)A' smk^ik) (27) 



?/'(even, even, even) 


= ^1 


^/'(odd, even, odd) 


= ip-j 


'?/'(even, odd, even) 


= V'e 


'(/'(even, even, odd) 


= 'ip5 


'(/'(odd, odd, even) 


= ^4 


^/'(odd, even, odd) ; 


= ^3 


'(/'(even, odd, odd) 


= ^2 


'(/'(odd, odd, odd) e 


= ^8 
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where 

= {ki : -7r/2 < A;^ < 7r/2} 
is the Brillouin zone of the (even,even,even) sublattice, 

A' = 

and 

i_/0 1\ 2_/cr^ 0\ 3_/ct3 q 









are a particular representation of the Dirac matrices. 
In this representation the mass operator is 



where 

B -- 

The Fermion spectrum is 



1 
-1 



3 

sin^ /cj + 



and only the region ~ is relevant to the continuum limit, 
normalized il^ik) so that 

{i^{x),i^\y)]^5{x-y) , ^ ^ik - I) 

If we define 

and the unitary matrix 



-a 



2\l + (3 1-/3 
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and 



iIj = Mtp' (37) 

with 

i^' = {^Pa,A) (38) 

the Hamiltonian is 

In the low momentum hmit, sin/cj ~ ki, with Fermion density 1/2 per site so 
that the Fermi level is at the intersection point of the positive and negative 
energy bands, we obtain 2 continuum Dirac Fermions. 

This describes two flavors of 4-component Dirac Fermions and the Dirac 
masses for each component given by the staggered charge density have op- 
posite signs. Thus the charge density breaks the discrete chiral symmetry. It 
also breaks a flavor symmetry which, in the absence of mass, mixes the two 
continuum Fermions. The continuum limit for arbitrary d is discussed in . 



3 Fermion Representation of SU(A^) Quan- 
tum Antiferromagnet 

In order to set the formalism which allows us to find the exact mapping be- 
tween strongly coupled lattice QED and SU(A^) quantum antiferromagnets, 
we review here the representation of an SU(A^) antiferromagnet in terms of 
fermion bilinears. The analysis holds for any space dimension D. 
The Hamiltonian for an U(A^) quantum antiferromagnet is 

^^AFM = ^ H Jab{x)Jba{y) (40) 
<a;,j/> 

where Jab{x), a,b = 1, ■ ■ ■ , A^l, obey current algebra relations associated with 
the Lie algebra of U(A^), 

[Jab{x), Jcd{y)] = {Jad{x)5bc " Jcb{x)5ad) - v) (41) 
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and where < x,y > denotes the hnk connecting sites x and ?/ on a bipartite 
lattice. For simphcity, we take the lattice to be cubic. Here, we have used a 
particular basis for the SU(A^) algebra which can be conveniently represented 
by Fermion bilinear operators, 

Jab{x)=rK^)i^''i^)-S"'/^ (42) 
with Fermion field operators obeying 

The representation of the algebra on each site x is fixed by specifying the 
Fermion number of the states, 

p{x)=Y.Jaa{x) (43) 
a 

For example, the Fermion vacuum state |0 > such that 

iP''{x)\0 >= , Va,x (44) 
is the singlet state, the states with m < N Fermions per site, 
nV^"'"^(a;)V^"'"^(a;)...V^"™^(x)|0 > 

X 

corresponds to p{x) = m — N/2 for all x and the irreducible representation 
with the Young Tableau with one column and m boxes. 

For each site x, p{x) is the generator of the U(l) subgroup of U(A^). Using 
a basis T' = (T*)*, ^ = 1, . . . , iV^ - 1, of the Lie algebra of SU(A^) in the 
fundamental representation normalized so that tr(T*T-') = T^fjT^^ = S^^ /2, 
and using 

Tabbed = ^^adSbc - ^SabScd (45) 

it is convenient to introduce 

r{x) = rHx)r,,r{x) (46) 

obeying current algebra of the Lie algebra of SU(A^), and to write the Hamil- 
tonian (^) as 

^AFM = T7 H Pix)piy) + HsviN) (47) 
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with 



Hsv(N) = g' E J'i^)-J'iy) (48) 

<x,y> 

is the Hamiltonian of an SU(A^) antiferromagnet. From this it is obvious 
that by fixing the -ffAFM is reduced to an SU(A^) antiferromagnet. 

For example, the famihar j = 1/2 SU(2) Heisenberg antiferromagnet is 
obtained from (0), = 2, by using the identity 

— — = -OadObc - -^OabOcd (.49) 

corresponding to P5| ) for N = 2. 

Generally, when N is even we will consider the representations where 
m = N/2, so that the Fermion occupation of each site is N/2 and p{x) = 0. 
When is odd we divide the lattice into two sublattices such that the 
nearest neighbors of all sites of one sublattice are in the other sublattice 
(when this is possible the lattice is said to be bipartite). When N is odd, 
the representation of SU(A^) has [N + l)/2 Fermions, i.e. p{x) = 1/2, on the 
sites of one sublattice and {N — l)/2 Fermions, i.e. p{x) = —1/2 on the sites 
of the other sublattice. 



4 QED on a Lattice 

In this Section we shall review the Hamiltonian formalism of Abelian gauge 
fields on a lattice. For the most part, this formalism can be found in some of 
the classic reviews of lattice gauge theory, for example. A novel feature 



of our approach is a careful treatment of the normal ordering of the charge 
operator and a discussion of the ensuing discrete symmetries. This normal 
ordering turns out to be important if the continuum limit has to have the 
correct behavior under C, P and T transformations. It will also be important 
in our later solution of the strong coupling limit. 

4.1 Hamiltonian and Gauge Constraint 

We shall discretize space as a cubic lattice and, in order to use the Hamilto- 
nian formalism; time is left continuous. Lattice gauge fields are introduced 
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through the hnk operators 



Ui{x) = e^^'(^) (50) 

which correspond to the hnks [x, i] of the lattice. Electric fields propagate 
on links of the lattice and the electric field operator Ei{x) associated with 
the link [x, i] is the canonical conjugate of the gauge field 

[Ai{x),E,{y)]=i6i,6{x-y) (51) 

The gauge field and electric field operators obey the relations 

A_i{x) = -Ai{x - i) , E_i{x) = -Ei{x - i) (52) 

The Hamiltonian of compact QED is 

H^c = T.jE!{x)+Y: ^,sm'{B[x,z,j]/2) 

[x,i] [x,i,j] 

+ E {tl^^r^x + Oe*^'(^V"(a;) + h.c.) (53) 

lx,i] 

where the second term contains a sum over plaquettes and the magnetic field 
is defined as the curvature of the gauge field at the plaquette [x,i,j], 

B[x,iJ] = A,{x) +Aj{x + l) + A_i{x + + A_j{x + j) 

= Ai{x) - Aj{x) + Aj{x + Ai{x+j) (54) 

As discussed in Section 2, the hopping parameters ^[i^^j] contains phases which 
produce a background magnetic fiux vr per plaquette. In the weak coupling 
continuum limit, the magnitude of tix,i] is one, |tp = = 1 in order that 

the speed of the free photon and free electron fields are equal, i.e. so that 
the low frequency dispersion relations for both the photon and electron have 
the same speed of light. However, in order to obtain a relativistic continuum 
limit in general it is necessary to make |t| a function of e^. We shall find that 
in the limit where is large, the speed of light is proportional to \t\/e and 
it is necessary that |t| ~ e as — > cxd. 

Equation ( p3D reduces to the standard Hamiltonian of QED in the naive 
weak coupling continuum limit. In the strong coupling limit compact QED 
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is confining and the phase transition associated with chiral symmetry 
breaking is generally of first order |26]. In (^31) we have introduced TVl flavors 



of lattice Fermions labelled by the index a = 1, . . . , Ni^. 

The Hamiltonian is supplemented with the constraint of gauge invar iance. 
The gauge transformations of the dynamical variables, 



Ei{x) 



A : 


Ai{x) 


A : 


Ei{x) 


A : 




A : 





are generated by the operator 

Gk = Y. (-VM^)E^{x) + A(x) [r\x)r{x) - NJ2 



(55) 



(56) 



The local generator of gauge transformations where A has compact support 
is 



A 



dA{x) 



(57) 



Both (|56D and (|57| ) commute with the Hamiltonians in (p3D. 

In (|56|) and (|57|) we have subtracted the constant from the charge 

density operator in order to make the gauge generator odd under the usual 
charge conjugation transformation 



A{x) 
E,{x) 



X) 



- Ai{x) 

- E,{x) 

(-l)Ef=i^'=^'^t(a,) 

(-l)Ef=i^fe^«(a;) 



(58) 



In fact, the Fermionic charge term in ( |5^ ) can be put in the manifestly odd 
form ^[^"•'^ {x) , if)"" {x)]. Of course, charge conjugation symmetry of the lattice 
theory is necessary to obtain charge conjugation of the continuum theory. 
We shall see later that, particularly at strong coupling, the subtraction term 
in (|56D plays an important role. It seems to have been ignored in previous 
literature (for example, see |2^, p5|). Its presence is particularly important 
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when A'l is odd since the charge operator has no zero eigenvalues in that 
case (the eigenvalues of il)°''^{x)ip°'{x) are integers). 

Chiral symmetry is related to translation invariance by one site. The 
Hamiltonian (|53D is invariant under the transformations 

Xj : Ai{x) Ai{x+j) 
Xj ■ Ei{x) Ei{x+j) 



Xj ■■ 



ip'^^x) (-l)S^=^+i'^'=^'*t(a; + j) (59) 



for j = 1, ■ ■ ■ , D. 

In the following we shall use the charge conjugation symmetry which is a 
combination of these two transformations: 



C ^ exi (60) 

This is necessary if the mass operators which we defined in the previous 
section are to be invariant under charge conjugation symmetry. Also, we 
shall see that the strong coupling ground state is invariant under C but not 
under either ^ or Xj by themselves. 

The dynamical problem of Hamiltonian lattice gauge theory is to find the 
eigenstates of the Hamiltonian operator ( ^31) and out of those eigenstates to 
find the ones which are gauge invariant, i.e. which obey the physical state 
condition (or, the "Gauss' law" constraint) 

g{x)\^p^ys. >= (61) 

In the case of compact QED, the periodicity of the gauge potential induces 
an additional symmetry which, being a large gauge symmetry, can be rep- 
resented projectively. In the Schrodinger picture we shall require that the 
quantum states which are functions of a configuration of the gauge field 
transform as 



\Ai{x) + 27m >= exp {ine[x, i]) \ Ai{x) > (62) 

There is a separate parameter 6'[x, for each link of the lattice. These pa- 
rameters originate in a way similar to the theta-angle in ordinary QCD. The 
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symmetry (^) together with the commutator (|5T| ) imply that the spectrum 
of the electric field operator has eigenvalues which are separated by integers 
and offset by 0: 

spectrum = {^^[x, i] + integers} (63) 

The Hamiltonian and gauge constraints can be obtained from the gauge 
invariant Lagrangian 

L = Y,r\x){idt - A^{x))r{x) + E e,{x)a,{x) 

^ [a;,«] 

+ J2 E^{x)V^Ao{x) + J2 Mx)Nl/2 - H (64) 

lx,i] X 

where the temporal component of the gauge field has been introduced to 
enforce the gauge constraint and the time derivative terms give the correct 
symplectic structure. In order to get Lorentz invariance of the Fermion spec- 
trum in the weak coupling (naive) continuum limit, we require half-filling of 
the Fermionic states, i.e. that the total charge defined by 

X 

has zero vacuum expectation value. 

Here we have considered massless QED. As well as the gauge invariance 
and charge conjugation invariance discussed above, the Hamiltonian is sym- 
metric under the discrete chiral transformations ( ^9|) which on the lattice 
corresponds to symmetry under translation by one site. In later Sections, we 
shall consider the possibility of spontaneous breaking of this symmetry. 



4.2 Gauge fixing and quantization 

We shall quantize the gauge fields in the Schrodinger picture. The quan- 
tum states are functions of the link operators which are represented by the 
classical variables Ai{x) and the electric field operators are derivatives 
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We shall also consider the usual Fock representation of the Fermion anticom- 
mutator. The empty vacuum is the cyclic vector 

^'^(a;)|0>=0 \/a,x (66) 

and Fermions occupying lattice sites are created by ip"'\x). 

In compact QED the spectrum of the gauge generator is discrete and 
a state which obeys the physical state condition can be normalized, thus 
implying that there is no need for additional gauge conditions. The basis 
wave-functions for compact QED (in the basis where the Fermions density 
and eletric field operators are diagonal) are ^'[n(a:)]$[A] with the Fermion 
states 

X a=l 

labeled by vectors n{x) = {ni{x), ■ ■ ■ ,niy^^{x)) with na{x) = or 1, and the 
photon states 

<^[A]=exp{iY,ei{x)A,{x)) (68) 

where the eigenvalues ei{x) of the electric field operator are in spectrum[i?i(a;)] 
(|63D. Furthermore the states of the photon field are normalized using the in- 
ner product 

< <^^[A],<^>,[A] >= n r^^^\[AmA] (69) 

[x,i] 

and the Fermion states have conventional inner product given by < 0|0 >= 1 
and the canonical anticommutator relations of the fermion field operators. 
The physical state condition ( pTf ) gives the additional restriction that 

V,e,(x) = -p{x) = -J2 EK(^) - 1/2) (70) 

X a 

where p{x) is the charge density (i.e. the eigenvalue of ip"''^ {x)'iIj°'{x) — Ni^/2). 
Pictorially, we can think of this as containing lines of electric flux joining sites 
whose charges are non-zero and also closed loops of electric flux. In the strong 
coupling limit, the Hamiltonian is diagonal in the basis (|68D. This gives a 
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natural explanation of confinement in compact QED in the strong coupling 
region. If we add a particle-antiparticle pair to a state in (^) it must be 
accompanied by at least a single line of electric flux. The energy of such 
a line of flux is proportional to its length. Therefore the electron-positron 
interaction grows linearly with distance and is confining. This is in contrast 
to the situation in non-compact QED where the electric flux is not quantized. 
In that particle-antiparticle pair can have many lines with arbitrarily 

small flux. The energy of the field is minimized by the usual Coulomb dipole 
configuration. In high enough dimensions this is not a confining interaction. 

5 Strong Coupling Expansion 

In compact QED the eigenstates of the electric field operator are normal- 
izable and can be used for the ground state. In this case we separate the 
Hamiltonian into three terms, 

lx,i] 

Hi = J2 t[-A {^^(^ + i)e''''^Mx) + h-c.) 

[x,i] 

^2= E ^sm'iB[x,t,j]/2) (71) 

In the strong coupling limit it is necessary to solve Gauss' law (^) for the 
electric fields and the charge distribution in such a way as to minimize Hq. 

When TVl is even, the charge operator has zero eigenvalues and Gauss' 
law has the solution where Ei{x) = 0, which is an obvious minimum of the 
Hq, and there are fermions on each site. The degeneracy of this state 

may be resolved , resulting in an effective Hamiltonian which describes the 
SU(iVL) antiferromagnet in the representation with Young Tableau having 
one column with boxes. This system has Neel order in D > 2 if A^l is 

small enough and the chiral symmetry of electrodynamics is broken, along 
with the SU(A^l) flavor symmetry [Q. 

In order to see this we follow [0 and resolve the degeneracy by diago- 
nalizing perturbations in the hopping parameter expansion. The flrst order 
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perturbations vanish. The first non-trivial order is second order, 

62 = - < g.s.\H. ^ H,\g.s. > (72) 

This matrix element can be evaluated by noting that Hi creates an eigenstate 
of Ho different from the ground states with additional energy 

AE = - + -p)Vi(x| ^|a;) = (73) 

2 2 ^ ^ ^ -V ■ V 

Diagonalizing the matrix of second order perturbations is equivalent to find- 
ing the spectrum of the effective Hamiltonian 

2\t\ 



Hes = ^5^^'^(x + z)V''(a;)r^(a;)V'"(a; + z) 

2\t " 



e2 ^ , 



e2 , .. 



Y.J^b{x)Jba{x + i) (74) 



where the operators Jab{x) = ip"''^ {x)4'^{x) — ^6ab, are the generators of the 
U(A^l) given in equation (^2|) and obeying the Lie algebra in equation (^T]). 
The constraint on the total occupation number of each site, 

Nl 

P{^) = Jaa{x) ~ 
a=l 

reduces to SU(A^l) and projects onto a fundamental representations of SU(A^l) 

i- 

It is straightforward to see that the higher orders in the hopping param- 
eter expansion also have higher orders of 1/e^. In fact, if we consider the 
following limit, 

— s> 00 , Itp 00 

|t|^/e^ = constant (75) 

all higher order perturbative contributions to both the wavefunction and the 
energy vanish. Thus, in this limit, QED is exactly equivalent to an SU(A^l) 
antiferromagnet. That (fTSD is the correct limit to take can be seen from the 
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fact that, if the antiferromagnet in ([7^ ) is in an ordered state, the speed of 
the spin-waves, which are the gapless low-energy excitations is proportional 
to \t\/e. They have linear dispersion relation uj{k) ~ \k\ and play the role of 
massless goldstone bosons for broken flavor symmetry. Their speed should be 
equal to the speed of light, which is one in our units. This implies that \t\/e 
should be adjusted so that the spin-wave spectrum is relativistic, uj{k) = \k\. 
Hence the limit in ([75|). 

When A^L = 2, this model is the quantum Heisenberg antiferromagnet in 
the j = 1/2 representation. It is known to have a Neel ordered ground state 



in D >3 p8[ and there is good numerical evidence that it has Neel order in 



D = 2. The antiferro magnetic order parameter is the mass operator 

S = ^(-l)I^f=i (76) 

which obtains a vacuum expectation value in the infinite volume limit. Thus, 
when Nl = 2 the strong coupling limit breaks chiral symmetry and generates 
electron mass. It is interesting that in this case there is an iso-vector con- 
densate. In the strong coupling limit this seems unavoidable. The only way 
to get an iso-scalar condensate is with a charge density wave. However such 
a configuration always has infinite coulomb energy compared to an electric 
charge neutral but isospin carrying condensate. 

The low energy excitations of this systems (with energies of order |tp/e^ 
are spin waves. All other excitations have energies which go to infinity in 
the limit (fTSl). The spin waves are the pions which are the scalar Gold- 
stone bosons arising from spontaneous breaking of the vector flavor symmetry 
SU(2)-.U(1). 

For large A'l there is some evidence that the SU(A^l) antiferromagnet in 
these specific representations has a disordered ground state [^. Particu- 
larly in 2+1-dimensions it is known that for infinite Nl the ground state is 
disordered 

When A^L is odd, since the charge density operator has no non-zero eigen- 
values, it is impossible to find a zero eigenstate of the Gauss' law constraint 
operator without some electric field. The problem which we must solve is to 
minimize the energy functional J2 E"^ subject to the constraint V ■ E = — p 
where at each site p has one of the values 

— —, — — + 1, . . . , — an odd mteger 
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The reader can easily convince herself that the charge distribution and elec- 
tric field which one obtains are [01, 



These configurations break chiral symmetry in that they are not invariant 
under the transformation xi in (H) but they are symmetric under charge 



conjugation C defined in (pO). 



The ground state energy per lattice site is 

\V\ 32D ^ ' 

Note that it is of order e^. This is in contrast to the ground state energy when 
the number of lattice Fermion fiavors is even, which is of order |tp/e^ ~ 1. 

The ground states that we have found are highly degenerate in that only 
the number of Fermions at each site is fixed. Their quantum state can still 
take up any orientation in the vector space which carries the representation 
of the fiavor SU(A^l) given by the Young Tableaux with one column and 
(A^L — l)/2 or (A^L + l)/2 rows on each of the two sublattices A and B in 
which we have divided the lattice. Namely : A is the set of all points where 

Xi is even and B the one where Y.i Xi is odd. The differing occupation 
numbers on sites on each sublattice yield different representations of SU{Nl). 

The degeneracy must be resolved by diagonalizing the perturbations, 
which are non-zero in second order and the problem is equivalent to di- 
agonalizing the antiferromagnet Hamiltonian ([7^). Also, in the limit (fTsP 
this correspondence is exact. 

The strong coupling ground states that we find when A^l is odd contains a 
charge density wave. The staggered charge density operator has expectation 
values 

-L < ^(-l)EL.^.^t(^)^(^) ±1 (78) 

This condensate is an isoscalar and we have shown that it must always occur 
in all dimensions. When A^l > 1 mass operators with certain generators of 
SU(A^l) could have expectation values if the ground state has antiferromag- 
netic order. However, unlike the case of even A^l, the antiferromagnetic order 
is not required in order to have chiral symmetry breaking. 
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The ground state we find breaks chiral symmetry. This is a true dynami- 
cal symmetry breaking since, in infinite volume, the ground states which are 
related by a chiral transformation are never mixed in any order of strong cou- 
pling perturbation theory. Furthermore, there are no local operators which 
couple them. 

We conclude that the strong coupling ground state breaks chiral symme- 
try for any odd A^l and in any dimensions. As in the case of even there 
is also the possibility (and for small A''^ the hkelyhood) that the SU(A^l) 
antiferromagnet we obtain here is in a Neel state and the fiavor symmetry is 



also broken 129 



Notice that in the ground state, the electric fields are not integers, but 
on each link, the spectrum of the electric field operator is 1/4D+ integers. 
The "theta angles" 1/4D survive all orders in strong coupling perturbation 
theory. 



6 Remarks 

In this paper we reviewed our approach to analyze the strong coupling 
limit of quantum electrodynamics using the Hamiltonian picture and a lattice 
regularization. We showed also the mechanism for chiral symmetry breaking 
in this lattice gauge theory. We chose to use staggered Fermions because 
they give the closest analog to interesting condensed matter physics systems. 

In dimensions, staggered Fermions give A^l species of 2- component 
Dirac Fermions. When A^l = 1 we obtain the Schwinger model with a lattice 
regularization. Also, in this case, we have found that the chiral symmetry is 
broken dynamically. Of course, due to the staggered Fermion regularization 
there is no continuous chiral symmetry, which is as it should be since it should 
be impossible to regularize the Schwinger model so that there is simultaneous 
continuous chiral and gauge symmetry. However, to match the solution of the 
continuum Schwinger model, the Fermion should obtain mass. This indeed 
happens in our strong coupling limit by spontaneous symmetry breaking. 
(Although we disagree with some aspects of the formalism, we agree with 



the results of reference |2^ on this point.) 

This result should not be confined to strong coupling, but should persist 
for all coupling, i.e. the critical coupling m D = 1 should be at = 0. We 
conjecture that this sort of symmetry breaking for small is a manifestation 
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of the Peierls instability — the tendency of a one dimensional Fermi gas 
to form a gap at the Fermi surface. This happens with any infinitesimal 
interaction. 

In fact, this must also happen for the case where A^l is even. Then, there 
cannot be any spin order in 1 dimension. However, anomalies break the 
isoscalar chiral symmetry in the continuum theory and should also do so here. 
This means that there should be a dynamical generation of charge density 
wave which would be driven by the Peierls insability. It also implies that for 
A^L = 2 for example, the ground state in the strong coupling limit would not 
be a Heisenberg antiferromagnet, but would be alternating empty site and 
site with two electrons in a spin singlet state. This state, even though it has 
large coulomb energy, avoids the infrared divergences of gapless Fermions. 

In higher dimensions, D > 2, one should explore the possibility of phase 
transitions between different symmetry breaking patterns for the SU(iVL) 
flavor symmetry as one varies A^l- There is already some work on this subject 
in the condensed matter physics literature on SU(A^l) antiferromagnets [p9| 



They analyze the SU(A^l) antiferromagnet which is similar to the strong 
coupling limit of an U(A^'c) gauge theory (see for details) which is in the 
representation corresponding to a rectangular Young Tableau with A^l rows 
and A^c columns. They work in the large A'c limit and show that there is a 
phase transition from the spin ordered Neel phase to a disordered phase when 
TVl ~ A'c. For QED A'c = 1 so their analysis is not accurate. For odd A^l the 
chiral symmetry is always broken and the question we are asking is whether 
the flavor symmetry is also broken. For even A^l possible phase transition is 
relevant to both chiral and flavor symmetry breaking. In this context it is 
interesting to ascertain if the order-disorder phase transition which occurs as 
one increases A^l in the SU(A^l) antiferromagnet is the same one that appears 
in the study of chiral symmetry breaking in 2+1-dimensional QED in the 



continuum [19, 20 1 where they find that chiral sjTumetry is broken only if the 



number of flavors is less than a certain critical value. These questions have 



been addressed in |jT^ 



Our study suggests that, besides the critical A^l, for a fixed A^l which 
is small enough, there should also be a critical couphng constant and, in 
fact, a critical line in the Ni^-e^ plane where there is a second order phase 
transition between a spin ordered chiral symmetry breaking phase and a 
disordered (and possibly chirally symmetric phase). We speculate that in 
3+1-dimensions a similar situation could occur. The latter phenomenon 
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has recently been discussed in the context of 2+1- dimensional gauge theory 
It can in principle shed light on the nature of the phase transition in 



the lattice theories. It would be interesting to analyze the critical behavior 
further, particularly in the lattice theory, since one could address the question 
of whether this phase transition is the same as the one seen in Josephson 
junction arrays. 

We have seen that the Coulomb gas representation of classical Josephson 
junction arrays is described at infinite coupling by the Nl = 1 Coulomb gas. 
Away from infinite coupling in the case of electrodynamics and away from 
the classical limit in the case of junctions arrays the two systems may still 
have many common features. 

In the quantum regime, Josephson junction arrays may be modeled by 
the quantum frustrated XY model which is equivalent to a model of lattice 
electrons in a magnetic field determined by the charge density [|T^. In a 
mean field theory, where the charge density is replaced by its average, this 
resembles the Azbel-Wannier-Hofstaeder problem for lattice electrons with 
nearest neighbor hopping with an external magnetic field |12|. When / = 
1/2 the low energy components of the electron spectrum are approximately 
described by a massless Dirac Hamiltonian. Thus, for / = 1/2, the model 
has relativistic lattice electrons and is gauge and parity invariant. When / is 
a rational number, the frustrated quantum XY model admits a flux phase as 



its mean field theory ground state |12|. It would be interesting to analyze the 



properties of this model beyond mean field theory. In particular to ascertain 
if - for / rational - the flux phase may be regarded as the exact ground state 
[0 as well as to investigate if, for / = 1/2, quantum fluctuations could 
reproduce the relevant operators of compact QED. 
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